Neel order, quantum spin liquids and quantum criticality in two dimensions 



Pouyan Ghaemi 1 and T. Senthil 1,2 

1 Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
2 Center for Condensed Matter Theory, Indian Institute of Science, Bangalore, India 560012 

This paper is concerned with the possibility of a direct second order transition out of a collinear 
Neel phase to a paramagnetic spin liquid in two dimensional quantum antiferromagnets. Contrary to 
conventional wisdom, we show that such second order quantum transitions can potentially occur to 
certain spin liquid states popular in theories of the cuprates. We provide a theory of this transition 
and study its universal properties in an e expansion. The existence of such a transition has a 
number of interesting implications for spin liquid based approaches to the underdoped cuprates. 
In particular it considerably clarifies existing ideas for incorporating antiferromagnetic long range 
order into such a spin liquid based approach. 
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I. INTRODUCTION 

In recent years a number of old ideas on theories of the 
cuprate materials have been clarified and rejuvenated. A 
variety of experiments now support the notion that the 
underdoped cuprates are usefully regarded as doped Mott 
insulators - in other words proximity to the Mott insula- 
tor strongly influences the properties of the underdoped 
materials^. The undoped Mott insulating parent mate- 
rials have long range antiferromagnetic order. However 
this order disappears rather quickly upon doping. Rather 
(at not too low temperatures) the pseudogap state that 
appears above the superconducting transition in the un- 
derdoped materials is probably best regarded as a doped 
version of a paramagnetic Mott insulator. Such an insu- 
lator has a 'built-in' spin (pseudo)gap. Vexing questions 
however remain on the precise theoretical connection be- 
tween such a point of view and the actual occurence of 
antiferromagnetic order in the undoped materials. 

A clue to resolving this dilemma is provided by neu- 
tron scattering experiments that reveal the existence 
of a sharp magnetic resonance at (jr, ir) in the doped 
superconductor 2 ' 3,4 . An appealing interpretation of this 
resonance is as a gapped version of the familiar magnon 
of a proximate antiferromagnetic stat o 5 i 6 i . Interestingly 
as the doping is reduced, the resonance frequency goes 
down proportionately to T^^. This suggests that if the 
doped state is to be viewed as a doped paramagnet, then 
the latter may at least be connected to the Neel state by 
a second order transition (see Fig^). Thus we are lead 
to search for quantum paramagnetic states of spin-1/2 
moments on a square lattice that are accessible from the 
collinear Neel state by a second order transition. 

It is important at this point to also review another 
old idea in cuprate theory. Early work 10,11 suggested 
that the undoped magnetically ordered Mott insulator is 
close to being disordered by quantum fluctuations into 
a paramagnetic featureless "spin liquid" phase. Such a 
spin liquid state was postulated to have neutral spin- 
1/2 spinon excitations and preserve all the symmetries 
of the underlying microscopic spin Hamiltonian (includ- 
ing spin rotation). Further it was argued that doping 




FIG. 1: Zero temperature phase diagram showing the route 
from Mott insulating antiferromagnet to d-wave superconduc- 
tor (dashed-dot line). Horizontal axis refers to doping and 
vertical axis refers to frustrating spin interactions that desta- 
bilize the Neel state. The thesis of the spin liquid approach is 
that the intermediate and long scale physics of the doped sys- 
tem may be fruitfully viewed as those of a doped spin liquid 
Mott insulator. Doping the spin-liquid phase naturally leads 
to d-wave superconducting state, (dashed lines) . 



a spin liquid could possibly lead to high temperature 
superconductivity 1011 . 

Despite its original appeal, this scenario was subse- 
quently questioned by a number of significant theoretical 
developments. Calculations in a controlled large- N ex- 
pansion of quantum Heisenberg spin models concluded 
that the natural result of destruction of collinear Neel 
order was not a featureless spin liquid but (for spin-1/2) 
a valence bond solid (VBS), which breaks various lat- 
tice symmetries^. (Here natural refers to the possibility 
that the quantum paramagnet in question is potentially 
separated from the Neel state by a second order transi- 
tion.) The VBS state also does not support fractional- 
ized spinon excitations. This was supported by a number 
of other indirect arguments - for instance by studies of 
quantum dimer models on the square lattice^. It was 
shown however that destruction of non-collinear Neel or- 
der could indeed lead to a fractionalized spin liquid state 
that preserves all lattice symmetries 1 ^. 

These calculations lead to the following folk wisdom 
(for a review see Ref. ITfil) : "In two spatial dimen- 
sions collinear ordered magnets naturally give way to 
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confined VBS paramagnets when disordered by quantum 
fluctuations while non-collinear magnets naturally lead 
to spin liquids" . As the magnetic ordering is undoubt- 
edly collinear in the cuprates this folk wisdom apparently 
spells doom for the view of the cuprates as doped spin 
liquid paramagnetic Mott states. 

The purpose of the present paper is to revisit these 
issues. We will first argue that, if at all a spin liquid 
based approach is to be pursued, experiments suggest a 
certain kind of paramagnetic spin liquid state as natural 
candidate 'parent' states of the doped cuprates. Next we 
argue that the existing theoretical work does not rule out 
a direct second order transition between the Neel state 
and this particular kind of spin liquid state. Finally we 
outline in some detail a theory for just such a direct sec- 
ond order transition. Thus our work calls into question 
the folklore described above and potentially frees the spin 
liquid based approach to the cuprates from one of its the- 
oretical criticisms. 

Based on these results we will develop a qualitative 
picture of the neutron resonance mode seen in experi- 
ments in the doped system and its relationship with other 
aspects of the observed spin physics. Our description 
will naturally unify two popular views of the resonance 
mode - one as a soft mode associated with antiferromag- 
netic long range order£&l, and the other as a spin ex- 
citon formed as a triplet particle-hole collective mode of 
fcrmionic quasiparticle9iS*i2^2ii2iSfi. 

We begin with experiments. It is by now quite clearly 
established that the cuprate superconductors are d-wave 
paired and furthermore have nodal BCS-like quasiparti- 
cles. The existence of the nodal quasiparticles is theoret- 
ically significant. Indeed the possibility of d-wave paired 
superconductors without nodal quasiparticles has been 
much emphasized by Kivelson and coworkers (for a re- 
view see Ref. l2lil . It is therefore of some interest to ask 
whether there exist paramagnetic Mott states that al- 
ready have gapless nodal excitations. Such a state then 
builds-in enough of the spin physics seen in the exper- 
iments at finite doping that it would be an attractive 
'parent' Mott insulator as a basis for a theory of the un- 
derdoped cuprates. 

Remarkably such states are known to exist as stable 
quantum phases^S of quantum antiferromagnets mag- 
nets on a two dimensional square lattice, at least within 
an appropriate large- N expansion. In this paper we will 
focus on one such state that has played a central role 
in some previous theoretical WO rkS4S5i2SiSLS2iS£ on the 
cuprate problem. This state - dubbed the ei-wave RVB 
or staggered flux (sF) spin liquid - is a quantum para- 
magnet that nevertheless has gapless spin carrying ex- 
citations. Recent theoretical worbSi has established the 
stability of such a state (in a suitable large-TV expan- 
sion) . A low energy description of the physics is usefully 
provided in terms of a theory of gapless nodal fermionic 
spinons coupled minimally to a fluctuating U(l) gauge 
field. Despite this however there really is no true quasi- 
particle description of the low energy spectrumS&SI. 



In this paper we are interested in exploring the quan- 
tum phase transition between this algebraic spin liquid 
and the collinear Neel state. Previous theoretical work 
providing the basis for the folk wisdom mentioned above 
do not constrain the nature of this transition. First the 
large- N calculations of Ref. were based on a bosonic 
representation of the spins. This representation is not 
well-suited to access the algebraic spin liquid phase. In- 
deed it is tailor-made to access either VBS phases or 
gapped spin liquids with bosonic spinons. Arguments for 
the presence of VBS order in the paramagnet based on 
quantum dimer models on the square latice 13 also do not 
help. Clearly the quantum dimer models are useful only 
for paramagnets with a full spin gap - and hence will not 
be able to access algebraic spin liquids. A more recent el- 
egant argument 31 attacks from the spin liquid side as fol- 
lows. First it supposes that to access Neel ordered states 
from spin liquids, the latter must have bosonic spinons. 
Then the Neel state is reached simply by condensing the 
bosons. Examining the dispersion relation of a bosonic 
spinon in a spin liquid state reveals that it generically has 
a minimum at an incommensurate wavevector. Condens- 
ing such a spinon then naturally leads to incommensurate 
spiral states and not to the simple collinear state. A loop- 
hole in this argument is the supposition that it is only a 
spin liquid state with bosonic spinons that can be prox- 
imate (i.e. separated by a second order transition) to a 
Neel state. Indeed we will show in this paper that the 
algebraic spin liquid state - which doesn't have bosonic 
spinons - can be connected to the Neel state by a second 
order transition. 

Our starting point is a mean field description of both 
the spin liquid and antiferromagnetic phases. This will 
be done in terms of a slave particle representation of the 
spin operator in terms of neutral fermionic spinons. At 
this mean field level the spin liquid state we study will 
have spinons that are paired into a d-wave state. The 
resulting spinon dispersion has four Fermi points in the 
Brillouin zone at which the spinons are gapless. Anti- 
ferromagnetism is obtained in this representation as a 
spin density wave transition of the spinons. This mag- 
netic ordering produces a gap for the spinons. Such a 
mean field description of the antiferromagnetic state was 
first proposed by Hsu2£. Similar ideas have subsequently 
been explor ed in a number of publications - see Refs. 
I3l3l35ll3l However it is important to realize that the 
resulting mean field state is not a conventional antiferro- 
magnet. Though the spinons have acquired a gap they 
still have not disappeared from the spectrum. Thus the 
mean field description is apparently that of a fractional- 
ized antiferromagnet. For the particular mean field state 
studied in this paper this problem is cured once fluctua- 
tions beyond the mean field are included. Indeed we will 
argue that these fluctuations confine the spinons in the 
magnetically ordered state. 

Does the slave particle mean field description of the 
antiferromagnet still have any physical meaning if the 
spinons are confined? An answer to this question is pro- 
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FIG. 2: Schematic phase diagram of quantum phase transi- 
tion between Neel and dRVB algebraic spin liquid phases. 



vided by considering the magnetic state close to the tran- 
sition to the spin liquid. The critical point we describe 
for this transition has some conceptual similarity with the 
deconfined quantum critical points^L^ studied recently. 
Specifically spinon variables will already prove useful in 
describing the critical point itself though they do not cor- 
respond to physical degrees of freedom deep in the Neel 
phase. Furthermore there are two diverging length/time 
scales - one of which diverges as a power of the other - 
as the transition is approached from the Neel side. The 
physics is that of the conventional Neel state only at the 
very longest scales. At length/time scales intermediate 
between the two diverging ones, the physics is correctly 
thought of in terms of the mean field spin density wave 
state formed out of fermionic Dirac-like spinons. This 
state (though ultimately unstable) describes an increas- 
ingly wider regime of intermediate scales close to the 
transition. Apart from the spin waves expected from 
broken spin symmetry it also contains gapped spin-1/2 
spinons coupled minimally to a gapless spin-0 gauge bo- 
son (a "photon"). Thus in this large intermediate scale 
regime the spinon mean field description of the antifer- 
romagnet is indeed the correct physical starting point. 
Note that accessing the true critical theory requires keep- 
ing fluctuations other than those of the Neel order pa- 
rameter. This should hardly be a surprise given the non- 
trivial nature of the paramagnetic phase. 

The rest of the paper is organized as follows. In Sec- 
tion [H] we quickly outline the mean field theory of the 
spin liquid and the transition to the Neel state. Next in 
Section mil we go beyond mean field, by including gauge 
fluctuation for spinons and dynamics for Neel field. We 
first briefly review the properties of the algebraic spin liq- 
uid phase with particular emphasis on precursor fluctu- 
ations of the magnetic ordering near the transition (Sec- 
tion [^^J. Then in subsection IIII Bl we discuss the mag- 
netic state and show that it gets smoothly connected to 
the conventional Neel state once monopole fluctuations 
are included. We very briefly discuss the connection to 
projected wavefunction descriptions of the magnetism in 
undoped spin models in subsection IIII CI Next in Sec- 
tion ]W\ we consider the phase transition. We show that 
within certain assumptions the transition has conceptu- 
ally similar structure to the deconfined critical points of 
Ref. |33- We then present results for both the ASL phase 



and the transition to the Neel state within an appropri- 
ate e expansion in Section ]V\ In Section fVTl we discuss 
implications of our results for theories of cuprates before 
concluding in Section IVIII Various appendices contain 
technical details of the calculations. 



II. MEAN FIELD THEORY 

Consider a generic SU(2) symmetric spin-1/2 model 
on a square lattice with predominantly antifcrromagnetic 
short ranged interactions: 

H = J & ■ S r > + ■■ ■ (1) 

<rr'> 

Here J > (antiferromagnetic exchange), and the ellip- 
sis represent frustrating interactions that can be used to 
tunc quantum phase transitions. We will require that the 
full Hamiltonian be invariant under SU(2) spin rotations, 
time-reversal, and the full space group of the square lat- 
tice. It is well known that the nearest neighbour model 
has a Neel ordered ground state. Various paramagnetic 
ground states can be accessed (in principle) by appropri- 
ate frustrating interactions. As explained in the intro- 
duction, here we will focus on a particular paramagnetic 
state that is known as the dRVB algebraic spin liquid 
(also often referred to as the staggered flux spin liquid). 
A mean field theory for this state has been described 
several times in the literature and is well-known2i. First 
the spin is formally rewritten as a bilinear of fermionic 
"spinon" operators 

Sr = l^fLvapfrp- ( 2 ) 

Here a = 1,2, corresponding to spin up/spin down 
fermions. This is an exact rewriting when combined 
with the local constraint f^f a = 1. In the mean 
field approximation the exact Hamiltonian is replaced 
by one quadratic in the spinon operators but with self- 
consistently determined parameters. For the dRVB state, 
the mean field Hamiltonian takes the form 

H sF = - J2 ((Xrr>+i&rr')fifr'+h.c) (3) 
<rr'> 

Here we take r to belong to one sublattice of the square 
lattice. So that r' belongs to the opposite sublattice. 
The constants A rr / = +A on horizontal bonds, and —A 
on vertical bonds while \rr> — t on all bonds. This de- 
scribes fermionic spin-1/2 spinons on the square lattice 
with complex hopping amplitudes such that there is a 
non-zero flux that is staggered from plaquette to pla- 
quette. Despite appearances, this saddle point possesses 
the full symmetry of the microscopic model including all 
lattice symmetries. (The apparent breaking of transla- 
tional symmetry is a gauge artifact). Recent work has 
clarified the nature of fluctuations about this mean field 
state. But in this present section, we will stick to the 
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FIG. 3: Two-site unit cells (indicated by dashed lines) used 
to diagonalize mean field Hamiltonian. 



mean field description and see how a transition to a Neel 
ordered state may be described. 

To access a Neel state we modify the dRVB mean 
field Hamiltonian by adding a nearest neighbour anti- 
ferromagnetic interaction between the spinons. Such an 
interaction will anyway be induced once fluctuations be- 
yond the mean field theory are considered. By including 
it explicitly, we can induce a spin density wave ordering 
of the fermionic spinons. We therefore consider 



{rr') 



'frfr' + T* r /f r /f r ) H S r .S r ' 



(4) 



(rr 1 ) 



Here T rr i = T = t + iA on bonds as shown in Fig. El and 
equals T* =f-iA on other bonds as also shown in Fig|3 
We now treat the ^ term in a mean field approximation. 

We look for a solution where < S r >= e r Nz is non-zero 
(In mean field theory, N = (e r S r ) is constant. So we can 
choose it's direction as z direction). Here e r = (— l) x+y 
is +1 on the A sublattice and —1 on the B sublattice. 
The mean field Hamiltonian reads 



H MF = ~ J2 (Tfrfr' + T*fl,f r ) - — 
(rr') 



2^ e rf r !y fr 



(5) 

The value of N is to be determined self-consistently. We 
can diagonalize this Hamiltonian using a two site unit cell 
as plotted in figure This gives the following equation 
for energy eigenvalues and eigenstates: 



A N2 1 e (k)-iA(k) 



e(k) + iA(k) 



g 2 



h 
h 



h 
h 



(6) 



Heres, f\ belongs to one sublattice and fi to the other 
one. We have decomposed T as t — iA. e(k) and A(k) 
are then defined as: 

e(k) = — 2t(cos(k y ) + cos(k x )) 

A(k) = 2A(cos(fc E ) - cos(k y )) 

This gives spectrum of the two bands: 
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A(k) 2 



(7) 



Now with this in hand we can derive the self-consistency 
equation for N: 



(- 



OH 



MF 



(8) 



Where L is the linear system size. Using the spectral 
function of the lower band (we consider T = 0) this gives: 



k 



(2N\ 2 



e(k) 2 + A(k) 2 



(9) 



This equation has a trivial solution N — 0, but it is 
obvious from Q that non-zero solution, if it exists, has 
lower energy. So the system has two phases, achieved by 
tuning the value of g. The critical value, g c is given by: 



(10) 



for g > g c there is no non-zero solution and so N = 0. 
But for g < g c we have N ^ 0. From © we can also 
derive the behavior of N at critical point within mean- 
field: 



N oc 



for g > g c 
9c - 9 for g < g c 



(11) 



Note that in the magnetic phase the non-zero N induces 
a gap to the spinons. 

To study low energy properties, we will set up a contin- 
uum effective theory. In the next section we will include 
fluctuations in this continuum field theory. In the spin 
liquid state, the spectrum consists of two Fermi points, 
located at k x = k y = % and k x = — k y = — ^ (at these 
points e(k) = A(k) = 0), in the reduced Brillouin zone. 
There are gapless spinon excitations near these nodes 
with a Dirac-like linear dispersion. A low energy de- 
scription of the spin liquid is then possible in terms of a 
continuum field theory of massless Dirac spinons (a brief 
review that helps fix notation is in appendix^) . To study 
the magnetic transition described above within this con- 
tinuum field theory, we need to introduce a 'mean field' 
that couples to the (tt, 7r) component of the physical spin 
density. The resulting action takes the form 

S m = J d 2 xdTi){-i^d^ + iX ^-N.^-)^ (12) 

In this representation, ip consists of four two-component 
Dirac fields. The four Dirac fields arise from the pres- 
ence of two physical spin species together with two pairs 
of nodes. The Pauli matrices o act on the spin index 
while the /2 are Pauli matrices acting on the node index. 
It is readily seem that the combination ^■ipii z <jip is pre- 
cisely the continuum form of the physical spin density 
near (ir,ir). In the mean field theory N is to be deter- 
mined self-consistently. The coupling A is proportional 
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to coupling i and from now on the momenta are consid- 
ered with respect to the nodes. As expected a non-zero 
N gaps out the Dirac spinons. This gap vanishes upon 
approaching the phase transition to the spin liquid. The 
inverse of this gap determines a diverging length scale - 
within the mean field theory this length scale describes 
the decay of the connected part of the spin correlations 
near (tt,tt). This may be seen by a direct calculation 
(described in appendix IO which gives: 

e^^(0)^(r)) c cx ^1(1 + Sij (13) 

So the connected correlation for r <C £ is a power-law 
decaying function with fourth power of r. For r 3> £ it is 
a exponential decaying function, with correlation length 
£ and a pre-factor which decays as the third power of r. 
The correlation length at critical point diverges as: 

|A-A C | \g-9c\ 

III. BEYOND MEAN FIELD THEORY 

In this Section we consider the effects of fluctuations 
beyond the mean field theory described above. In the 
spin liquid phase far from the magnetic transition the 
crucial fluctuations are those associated with the phase of 
the spinon hopping parameter. These are to be thought 
of as gauge fluctuations associated with a (compact) U(l) 
gauge field that is coupled minimally to the spinons. Re- 
cent work has shown that the dRVB spin liquid is sta- 
ble to such gauge fluctuations 23 (at least within a sys- 
tematic large-iV expansion where N is the number of 
Dirac spinons). Through out this paper we will assume 
that this stability persists to the physically relevant case 
N = 4. The low energy theory of the resulting phase is 
described by massless QED in three space-time dimen- 
sions: 

S = J d 2 xdr {-#7^ + i e o M )]V + (e^ K d u a K ) 2 } 

(15) 

Here a M is a fluctuating U(l) gauge field which may be 
taken to be non-compact at the low energies. This the- 
ory flows to a conformally invariant fixed point. Various 
physical quantities have non-trivial power law correla- 
tions at the resultant spin liquid fixed points^. In par- 
ticular the (tt, 7r) spin correlator decays as a power law: 

e^(S(0).S(r) ~ -L (16) 

The exponent A is not known - a rough estimate from 
projected wavefunction o 40 ! 41 gives A w 0.75. The dy- 
namical spin correlations at (71", tt) in the scaling limit 
follow straighforwardly from the relativistic invariance of 



the field theory above. For the full zero temperature dy- 
namical spin susceptibility, we have 

xH,L(q,»)~-LF(-) (i7) 

w z " \vqj 

Here q is the deviation of the wavevector from Q = (tt, tt), 
F is a universal scaling function, and v is a non-universal 
spinon velocity associated with the nodal Dirac disper- 
sion. The exponent n is the anamolous dimension of the 
staggered spin and is related to A through 2A = 1 + 77. 
Due to the power law spin correlations, this spin liq- 
uid phase has been dubbed as "algebraic spin liquid" 
(ASL) 3 — . Note that the spinons are not good quasiparti- 
cles at low energies in the spin liquid phase. Indeed there 
presumably is no quasiparticle description of the spec- 
trum (rather like at interacting quantum critical points) . 
Nevertheless the field theory above in terms of spinons 
provides a useful description of the system. 

A remarkable feature of the dRVB algebraic spin liq- 
uid phase is the emergence of a huge global symmetry 
group characterizing the low energy fixed point. The low 
energy theory has an SU (4) symmetry corresponding to 
free unitary rotations between the four Dirac species. In 
addition the irrelevance of space-time monopoles at low 
energies implies a non-trivial global U(X) symmetry asso- 
ciated physically with the conservation of internal mag- 
netic flux. Ref. |3(J studied a number of consequences of 
the SU(A) symmetry. In particular it showed that several 
other competing order parameters had the same power 
law correlators as the Neel vector - these include the or- 
der parameter associated with the columnar/plaquette 
VBS orders. 

Near the transition to the antiferromagnetic state, we 
must treat the Neel field introduced in the previous sec- 
tion as a fluctuating vector N. Further upon integrating 
out high energy spinons (i.e. ones far away from the 
nodes), this N field will develop some dynamics of its 
own. The resulting action takes the form 

S = J d 2 xdr {-#7 M (d M + ie a„)]if> + iA$(/x* N.^)ip+ 

(18) 

In writing this action we have ignored anisotropics in 
the spinon velocities at the Dirac node and any difference 
between the velocities of spinon and N fields. Later we 
will show that all of these velocity anisotropies are irrel- 
evant (if small) at the critical fixed point between the 
dRVB ASL and Neel states (see EH . 

In the presence of the coupling to the N field, the ac- 
tion no longer has full SU (4) global symmetry. An SU (2) 
subgroup - corresponding to physical spin rotations - is 
still obviously a symmetry. This involves an SU(2) spin 
rotation of the i/j field together with an 0(3) rotation of 
the N vector. In addition, the global transformation 

^-fe^V (19) 
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with 9 a constant is also a symmetry. This is a U(l) sub- 
group of the full SU (4) symmetry. Thus the action has a 
global SU(2) x U(l) symmetry apart from the Un ux (l) 
associated with the gauge flux conservation. The extra 
U(l) symmetry that survives from the full SU(4) has the 
consequence that fermion bilinears such as ip(i x ?p,?pfj,y'ip 
can be freely rotated into another. In the original spin 
model, these operators transform identically to the VBS 
order parameter. The U(l) symmetry then implies that 
the columnar and plaquette order parameters can be ro- 
tated into one another, and hence will have identical cor- 
relations. 

Let us now study some general aspects of the two 
phases, near the phase transition. 

A. Precursor fluctuations in the spin liquid 

We will first consider the precursor fluctuations of the 
magnetic ordering in the spin liquid side. 

First consider the limit A = 0. Then the N field de- 
couples from the spinon-gauge sector. It is instructive 
to think about the spectral function for the (tt, tt) spin 
correlations in the spin liquid phase in this limit. It is 
simply a sum of two pieces as shown in Fig. 0]- a diverg- 
ing power law coming from the ASL, and a sharp delta 
function peak coming from the N fluctuations. Now con- 
sider turning on a small non-zero A. The low frequency 
divergence of the spin susceptibility will be unaffected by 
the coupling to the fluctuating N field (this follows from 
the assumed stability of the ASL fixed point). However 
the delta function peak coming from the N field will now 
be broadened due to decay into two spinons. The broad- 
ening may be described within a simple RPA approxima- 
tion (for details see appendix ITj|) . 
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FIG. 4: Dynamical spin susceptibility at (tt,tt) in the spin 
liquid phase in A=0 limit 

In Fig we plot the dynamical spin susceptibility at 
(tt, tt) as the phase transition is approached (by decreas- 
ing r) . Note that as expected, the peak coming from the 
N fluctuations "softens" on approaching the transition. 



B. Magnetic state 

When the parameter r is sufficiently negative the N 
field will condense leading to magnetic long range order. 
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(c) 

FIG. 5: Dynamical spin susceptibility at (tt, tt) after turning 
on a nonzero value for A. The plot sh ows t he change upon ap- 
proaching the transition from |5(a'J| to |5(cj| Note the softening 
of the N peak, as the transition occurs. 



In the continuum field theory this may be viewed as a 
spin density wave state that arises out of the dRVB ASL. 
In this subsection we will argue that contrary to naive 
expectations, it is in fact a different state from a conven- 
tional Neel state - rather it is a fractionalized antiferro- 
magnet in the same spirit as that studied in Refs. 12^1331 
This apparent problem will be cured once we include the 
effects of monopole fluctuations (ignored so far). We will 
show that this fractionalized antiferromagnet evolves at 
long length/time scales into the conventional Neel anti- 
ferromagnet. 

Consider first the description of the Neel ordered state 
within the mean field theory developed in Section [H] 
The mean field spectrum consists of gapped spin-1/2 
spinons 51 . It is important to realize that the spinons 
are merely gapped - they however have not disappeared 
from the spectrum. Now consider including fluctuations 
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as in Section ITTT1 The important fluctuations are those 
associated with slow rotation of the direction of the Neel 
order parameter N (spin waves) and those asociated with 
the phase of the fermion hopping T (gauge fluctuations). 
These are both conveniently discussed within the contin- 
uum theory in Eqn. El that obtains close to the critical 
point. Integrating out the gapped spinons, we may ob- 
tain an effective action for the spin waves and the gauge 
fluctuations. To quadratic order in both the transverse 
component of the Neel vector N± and the gauge field a, 
we get 




where the ellipses refer to higher order terms that are 
unimportant at low energies. The first term describes 
the expected gapless spin wave excitations. The second 
term describes a gapless linear dispersing "photon" as- 
sociated with the gauge fluctuations. This extra gapless 
mode provides a sharp low energy distinction between 
this Neel state (as described so far) and the conventional 
one. The gapless photon mode is minimally coupled to 
the gapped spinons - the presence of gapped spinons 
serves as another distinction with the conventional Neel 
state. Thus the antiferromagnet state is to be character- 
ized as a "fractionalized antiferromagnet" with a U(l) 
gauge structure. Following the notation of Ref. 133. we 
will dub it 17(1) AF*. 

Let us now include monopole fluctuations. In this 
magneticaly ordered phase, the low energy gauge ac- 
tion is that of free Maxwell theory in 2 + 1 dimensions. 
Then standard arguments show that the monopoles are 
strongly relevant. Thus the i7(l) AF* state (in two di- 
mensions) is ultimately unstable to monopole prolifera- 
tion. The result is to gap out the photon mode and cause 
confinement of all objects that carry non-zero gauge 
charge. In particular it implies that the spinons (which 
survived as gapped excitations when monopoles were ig- 
nored) will now be confined and disappear from the spec- 
trum. The resulting state is thus simply smoothly con- 
nected to the conventional Neel state. Thus including 
monopole fluctuations cause an instability of the uncon- 
ventional U(l) AF* state toward the conventional Neel 
state. 



C. Projected wavefunctions 

Before continuing we digress briefly to make contact 
with the large body of work on Gutzwiller projected su- 
perconducting wavefunctions (see Ref. and references 
therein), as a route to implementing RVB ideas. Of inter- 
est to us, will be studies on projected d-wave BCS states 
and their variants. In the slave particle description, a 
useful guess for a prototypical wavefunction for a state 
is obtained by taking the mean field state and projecting 
it onto the space of physical states. At half-filling this is 
equivalent to doing a Gutzwiller projection on the mean 



field state. According to this prescription, a guess for 
the wavefunction of the dRVB algebraic spin liquid will 
simply be 

\dRVB >= P G \dBCS > (21) 

where \dBCS > is the mean field ground state of a d-wave 
superconductor at half-filling with just nearest neighbour 
hopping and pairing on the square lattice. Correspond- 
ingly, a guess for the wavefunction of the magnetic state, 
as we have obtained it, would simply be 

| AF >= P G \dBCS + SDW > (22) 

The preprojected state on the right simply has spin den- 
sity wave order at (it, w) coexisting with the d-wave su- 
perconductivity. Such wavefunctions have been studied 
numerically 42 and are known to have excellent energy for 
the nearest neighbour Heisenberg model. From our con- 
siderations in the previous section, we would expect that 
this wavefunction is a prototype for a confined antiferro- 
magnet with no finite energy spinons. Some support for 
this expectation comes from the work of Ref. |43j which 
studied the properties of a single hole in that state. The 
quasiparticle residue was found to be non-zero consistent 
with that expected in a confined antiferromagnet. 

IV. PHASE TRANSITION: GENERALITIES 

Let us now consider the phase transition between the 
dRVB ASL and the Neel state. In the limit A = 0, the N 
vector fluctuations are decoupled from the ASL and the 
magnetic transition is simply in the universality class of 
the usual 0(3) fixed point in D = 3 space-time dimen- 
sions. Note that in this limit mean field theory predicts 
that the Neel order parameter vanishes with exponent 
(3 = 1/2 on approaching the transition. What is the ef- 
fect of turning on a weak A at this decoupled transition? 
First note that in the mean field theory of Section ITT1 we 
found that the Neel order parameter vanished with expo- 
nent (3 = 1 clearly different from the A = limit. Thus a 
non-zero A already changes the answers within mean field 
theory. More generally the effects of a weak A may be 
assessed by considering the renormalization group flow of 
A at the decoupled fixed point. We have 

d\ , . < , 

— = (D - A N - A)A (23) 

where D = 3 is the space-time dimension, A at is the 
scaling dimension of the N field at the D = 3 0(3) fixed 
point, and A is the scaling dimension of the spin oper- 
ator near (tt,tt) at the dRVB ASL fixed point. Here I 
is the usual logarithmic renormalization scale. We have 
A N = 1£24L a i, and A = ^±2. Thus D — A N -Aw 
2 — With the rough estimate rj w 0.5, we find that A 
is strongly relevant at the decoupled fixed point. Thus 
the true critical behavior will involve strong coupling be- 
tween the N field and the spinons of the ASL. In Section 
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fVl we will study this critical behavior in a controlled 
3 — e dimension. A very similar field theory where a fluc- 
tuating 0(3) vector field was coupled to massless Dirac 
fcrmions was studied many years ago by Balents et al 4 ^ 
in a different physical context. The main difference be- 
tween the theory of Balents et al and the action in Eqn. 
I18lis the presence of the gapless gauge fields in the latter. 
We will see that this modifies the universality class of the 
transition from that in Ref. 

What about monopole fluctuations at this critical 
point? Let us first review the situation in the para- 
magnetic algebraic spin liquid state. Recent work has 
argued 23 that when the number of Dirac species N is suf- 
ficiently large (i.e. bigger than some N c ) the monopoles 
are formally irrelevant at the ASL fixed pointS. In the 
large- N expansion, the monopole scaling dimension will 
be o(N) both at the ASL fixed point and at the crit- 
ical fixed point - thus at least for large enough N the 
monopoles are irrelevant at the critical fixed point as 
well. In this paper we will make the crucial assumption 
that this irrelevance continues to hold at N = 4, i.e. for 
SU(2) spin models (see figure El- 



Critical U(l) AF* 
\—^ 



AF 



s 



o 
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U(l) AF* gc ASL 

{Critical Coupling) 

FIG. 6: Renormalization flow diagram near the critical fixed 
point. The vertical axis is the monopole fugacity; the hori- 
zontal axis is a coupling g which describes the strength of the 
short range part of the spinon interaction. 

With this assumption the monopole fugacity is irrel- 
evant at the critical fixed point (and the paramagnetic 
ASL fixed point) but relevant at the ordered fixed point 
of the continuum field theory in Eqn. 1181 In renormaliza- 
tion group language, the monopole fugacity is a danger- 
ously irrelevant coupling. The length scale £ m at which 
the photon gets gapped (which may loosely be dubbed 
the "confinement scale" ) in the magnetic side may be es- 
timated as follows. Let the monopole scaling dimension 
at the critical point be A m > 3. Upon scaling out of 
the critical region in the ordered state to the correlation 
length scale £, the monopole fugacity will renormalize to 
~ £3-A m _ j-j. j g beyond this scale £ that the action 
in Eqn. [2U| starts applying. In the free Maxwell theory 
that obtains beyond £, the monopole fugacity grows. A 
standard matching argument now gives £ m ~ £(A m -i)/2 
Thus £ m diverges faster than £. The physics on scales 
£ <C L <C £ m is that of the fractionalized antiferromag- 
net 17(1) AF* . It is only at the longest scales L ^s> £ m 
that the conventional Neel behavior is obtained (fig. 0). 
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FIG. 7: Crossover length scales in the magnetic state close 
to the transition to the spin liquid. The shorter length scale 
£ describes the crossover from the critical state to the frac- 
tionalized antiferromagnet. The longer scale £ m is where this 
exotic antiferromagnet crosses over to the conventional Neel 
state through confinement. Both scales diverge near the crit- 
ical point but £ m diverges faster than f. 



V. e EXPANSION FOR CRITICAL 
PROPERTIES 

In this section we will show how the structure of the 
critical fixed point can be studied in a formal expansion 
near three space dimensions. Some care is necessary in 
dealing with the Dirac matrices in arbitrary dimension. 
But here following Ref. 0, we will sidestep this issue. We 
will perform calculations in a perturbative expansion of 
the coupling constants in d = 2, take the traces over the 
Dirac matrices and then finally, in evaluating momentum 
integrals, set d = 3 — e. As a warm-up we first describe 
the ASL fixed point within this e expansion. 



A. ASL in the e expansion 

Before starting the e expansion studies for the full the- 
ory, we study ASL phase with this approach. To do that 
we consider action (|18|l and assume N = 0. Then it 
reduces to the usual QEDj, action. The flow equation 
for this theory is well knowni 5 . (and also calculated in 
appendix Q : 



ee 



16 e 4 

y (i^f 



(24) 



Here e = 3 — d. This flow equation indicates the presence 
of nontrivial fixed point of order e at: 



*2 o 2 



(25) 



This is the pure ASL fixed point. The microscopic deriva- 
tion of the continuum field theory for the spin liquid al- 
lows for a velocity anisotropy between the two spinon 
nodes. This anisotropy was found to be irrelevant in the 
large- N limit of QED^^^. Here we examine its fate 
within the e expansion. Direct calculation fappcndix IC 5|) 



5 e, where S mea- 
= corresponds to 



shows that (3s = — ^ gpfp e* 2 = — | 
sures the velocity anisotropy (i.e. S 
isotropic QED3). So that it is irrelevant at o(e). Com- 
bined with the large- N result, this is strong evidence for 
its irrelevance at the physically relevant ASL fixed point 
for N = 4 in two space dimensions. As noted in previ- 
ous papers^ this irrelevance implies that the N = 4 ASL 
fixed point has global SU (4) symmetry corresponding to 
free unitary rotation between the four Dirac species. 
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Finally we can examine the scaling of gauge neutral 
fermion bilinears (such as the Neel vector) in the e- 
expansion. This is conveniently done by adding a source 
term that couples to such a bilinear, and calculating the 
one loop correction to the corresponding vertex (see ap- 
pendixEJ- We find that the (it, ir) component of the spin 
has scaling dimension A = 3 — 1.94e. Setting e = 1 gives 
the estimate A w 1.06. This implies extremely slow de- 
cay of the corresponding correlator. This estimate may 
be compared with that from the 1/N expansion directly 
in d — 2 which gives A as 1.54. Thus both expansions 
give slow decay for the Neel correlations that are strongly 
enhanced compared to the mean field results. 



B. Critical fixed point 

In this section we study the critical point within the 
e expansion by including the fluctuating N field. ^From 
the action IjlSjl we have three types of vertices and eleven 
different one loop diagrams. At one loop level, as derived 
in appendix[21 we get the following set of flow equations: 



equations admit another fixed point located at: 



Pu 
Pr 



16 



ee 



eX 2 



eu 



3 (4tt) 2 



10 

" (I^j 2 
11 u 2 



A 4 + 



10 

(4^)2 

\ 2 u 



e 2 X 2 



3 (4tt) 2 
u 



16 



(4tt) 2 
A 2 



96 



(4tt) 2 



3 (4tt) 2 (4tt) : 



(26) 
(27) 
(28) 
(29) 



A* = 



37r 2 e 
23^ 2 



12V 



(-36 + Vl2934)e 



(33) 
(34) 

(35) 



It is readily checked that the e 2 = fixed point is unsta- 
ble towards this one. Thus the presence of the gauge field 
has changed the universality class of the transition. Here 
we should also note that the at the one loop level, calcula- 
tion in appendix IC 5l shows that at this fixed point, the ve- 
locity anisotropy is irrelevant at o(e). So this fixed point 
is also stable against small velocity anisotropy. Now us- 
ing the flow equation for r (Eqn. I29f) . it is easy to extract 
exponent v for this fixed point: 



- = 2 - 4.07e 
v 



(36) 



Note that simply setting e = 1 gives an unphysical an- 
swer. This is a signal that the leading order e expansion 
is not quantitatively very accurate in estimating scaling 
dimensions in two space dimensions. Despite this the e 
expansion is useful to describe the structure of the fixed 
points and the trends of the various exponents. 

It is very interesting to ask about the behavior of the 
staggered spin correlations (i.e. near (n, 7r)) at this criti- 
cal point. Naively there are two different physical opera- 
tors that have the same symmetries as the staggered spin: 
the vector N and the fermion bilinear 2Vj| = ip(A z aip. 
Thus in writing down an expression for the staggered 
spin in terms of the fields of the continuum theory, we 
must include both contributions: 



Note that flow equation for electric charge is the same 
as usual quantum electrodynamics. In fact, gauge in- 
variance dictates this iormM^ (to keep the form e a M 
invariant under RG flow we need to have Z e / 1 \fZ~ a = 1)- 
From these we can get the following fixed points: 



JQ-S 



S(x) ~ Cl N + c 2 N z A + 



(37) 



,*2 



i *2 





8^ 
5 6 

384tt 2 
55 



(30) 
(31) 
(32) 



This fixed point describes the transition in the absence 
of the gauge field and was first discussed by Balents et 
al-Ll. Our calculations at this fixed point matches this 
previous work which therefore provides a useful check. 
As another check at this fixed point, if we consider our 
theory, with the three component N field replaced by a 
scalar field </>, it represent a Yukawa like theory, which has 
been studied in Ref. Our results then can be partially 
checked against these previous calculations. The full flow 



The ellipses refer to other operators with larger scaling di- 
mension that also have the same symmetries as the stag- 
gered spin. The coefficients are non-universal. Now 
consider the scaling of the staggered spin. If A = 0, then 
the fermion bilinear and N scale independently. Near 
d = 3, and with the available estimate of the scaling di- 
mension of the fermion bilinear at the ASL, it is readily 
checked that N has the lower scaling dimension. Hence 
the long distance decay of the staggered spin correlations 
will be determined by N in the A = limit. What hap- 
pens when A is non-zero as at the non-trivial fixed point 
above in the e expansion? It is expected that the true 
scaling fields will be some linear combinations of N and 



iVj| which will have the form 



AiN ■ 



A 2 NX 

B 2 m 



(38) 
(39) 



These fields will have scaling dimension Ai^ with (by 
definition) Ai < A2 so that the long distance decay will 
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be dominated by 4>\. The coefficients A\$ and B1.2 will 
be determined by the fixed point theory. At o(e) we ex- 
pect Ax ~ o(l),A 2 ~ o(e),5i - o(e), J B 2 - o(l). Thus 
to obtain the scaling dimensions Aj. 2 to o(e) we can ig- 
nore the 'mixing' terms A%,Bi, and simply calculate the 
anamolous dimension of N and N%. 

The exponent 77 for iV field (which determines the scal- 
ing dimension) is easily calculated from the field renor- 
malization coefficient Zn' 



7] is then given by coefficient of - in Z^: 

i] = 2.3e (41) 

so that Ai = 1 + 0.65e. The dimension A2 is readily 
calculated as in our discussion of the ASL above. We 
find A 2 = 3- 1.65e. 

It is straightforward to determine the scaling dimen- 
sion of all the fermion bilinears related to N A by SU (4) 
rotations. These are listed in Table [U the correspond- 
ing Feynman diagrams are in Appendix |0 The absence 
of ST/ (4) symmetry at the critical fixed point implies 
that these bilinears mostly all have different scaling di- 
mensions. Some weak constraints follows from the U(l) 
subgroup of the SU(4) that remains unbroken. For in- 
stance it implies that N^,N^ have the same scaling di- 
mension. As emphasized before physically this implies 
identical scaling of the plaquette and columnar VBS or- 
der parameters at this critical point. 



Field 


Spin Model 


Scaling 


theory 




dimension 


N%,N y A 


(-l) r - +1 & X S r + y , (-1) P »& X S r+3 


3-0.5e 


NX 


(-l) r * +r »5V 


3-1.65e 


N B 


(-l) r - +r * [(& + &)(&•&) 






+02 + Si)0i ■ S3)] 


3-1.65e 


Nh,N y c 


f 1) V Sr ' Sr-\-y j ( l) Sr ' Sr-\-x 


3-2.8e 


N z c 


[Si ■ 2 x Si) - S 2 ■ (S3 x Si) 






+ S3-(SiXS 2 )-S4-(Si xS 3 )] 


3+0.65e 


M 


[Si ■ 02 X Si) + S 2 ■ {S3 X Si) 






+ S 3 ■ 0i x S 2 ) + S4 ■ 0i xS 3 j\ 


3+0.65e 



-» 2 

We see that Na corresponds to Neel vector and iN^+N]! 
corresponds to VBS order parameter. 

C. Discussion 

Now let us examine the trends shown by the exponents 
calculated in table [I] Note that the scaling dimension 
of these gauge invariant bilinears are the same in ASL 
phase due to the SC/(4) symmetry. So for all of them 
we have A = 3 — 1.94e. Remarkably the scaling dimen- 
sion of the VBS order parameter (Nq, Nq) is smaller at 
the critical point than it is in the ASL phase. Thus the 
VBS fluctuations are enhanced by the critical N vector 
fluctuations. All other fermion bilinears decay faster at 
the critical point as their scaling dimension is increased. 
This includes the vector Na which is the contribution 
from the gapless fermions to the Neel vector. It is at 
present not clear whether in d = 2 the susceptibilities of 
these other operators (such as N^, etc) will diverge at 
the critical point (though they apparently do in the ASL 
phase). In contrast the VBS susceptibility will presum- 
ably diverge at the critical point. The divergence will be 
faster at the critical point (as say a function of temper- 
ature) than in the ASL phase. On the magnetic side the 
VBS susceptibility will of course be finite in the ground 
state. However it will diverge as the transition is ap- 
proached and will thus be large if the antiferromagnet is 
to be regarded as being close to this critical point. Thus 
a qualitative conclusion from our calculations is that in 
the limit that the Neel state can be usefully regarded as 
being born out of the dRVB spin liquid, it will also have 
enhanced VBS susceptibility. 

The diverging VBS susceptibility also provides an in- 
teresting way to define the corelation length £ in terms of 
directly measurable quantities. Consider the VBS corc- 
lations in the magnetic side as a function of length scale. 
At scales smaller than £ they will decay as a power law. 
However at scales larger a length set by £, they will de- 
cay exponentially. Thus £ may be usefully defined as the 
correlation length for VBS fluctuations in the ordered 
antiferromagnet . 

VI. IMPLICATIONS FOR CUPRATE THEORY 



TABLE I: List of observable in the spin model that are 
symmetry-equivalent to the N a and M fermion bilinears. For 
some of these we label the sites around the plaquette with 



lower-left corner at 

Si = Sr, S 2 = S r + x> 



f by the numbers 1, 
S3 = Sr+x+y and £4 



4. Precisely, 



r + y ■ 



Here we have defined the observable as: 



Na* = 


—iip/j, l aip 


(42) 


N B = 


—iipaip 


(43) 


N c l = 


—iip/i 1 -^ 


(44) 


M = 


—ii^ip 


(45) 



We now discuss some of the implications of our results 
for theories of the cuprates. As we emphasized in the in- 
troduction, a second order transition between a collinear 
Neel state and a gapless spin liquid is attractive for a 
number of reasons. Here we explore this in some greater 
detail. Our thinking on the cuprates is guided by Fig 
n We suppose that increasing magnetic frustration (the 
parameter g) at zero doping can induce a transition out 
of the collinear Neel state to a spin liquid. Theoretically 
the spin liquid is expected to evolve rather naturally into 
a superconductor when it is doped. The real material 
starts off in the antiferromagnetic state at zero doping. 
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The idea is that doping (apart from introducing holes) 
also has the effect of increasing g. Then we can hope that 
the intermediate and long scale physics of the resulting 
doped superconductor may be fruitfully described as a 
doped spin liquid. This is the rationale behind the spin 
liquid based approach to the cuprates. 

With this point of view in mind let us consider the 
effects of doping the dRVB algebraic spin liquid phase 
discussed in this paper. Previous papers (for a review 
see Ref. IHof) have shown how a d-wave superconductor 
with gapless nodal quasiparticles emerges quite naturally 
upon doping this spin liquid. Now consider reducing the 
doping in the real material. According to Fig Q] this also 
has the effect of reducing the magnetic frustration g. This 
pushes the "parent" spin liquid state closer to the tran- 
sition to antifcrromagnetism. The magnetic response of 
the parent spin liquid at wave vector (tt, tt) then evolves 
in the manner shown in Fig [5] Note in particular that 
if the transition to the magnetic state is second order 
then the "resonance" due to the N fluctuations softens. 
How does this impact the magnetic response in the doped 
superconductor? 

The doping of the spin liquid is incorporated theoreti- 
cally by the introduction of two species of charged spin- 
less bosons bi and 62- These bosons also carry gauge 
charges +1 and —1 respectively. Superconductivity is 
achieved when both bi and b 2 condense with equal am- 
plitude. This route from the dRVB spin liquid to the 
dSC has two important features. First the gauge charge 
carried by the bosons implies that the gapless gauge fluc- 
tuations of the spin liquid are quenched in the supercon- 
ducting state (by the Anderson Higgs mechanism). The 
spinons evolve naturally into the fermionic quasiparti- 
cle excitations of the dSC. The nodal structure of the 
spinons is retained - however coupling between the b and 
/ fields moves the nodes of the quasiparticles away from 
(7r/2, 7r/2) by an amount proportional to the doping x. 

For the magnetic response this has some crucial impli- 
cations. First when compared with Fig [5] the diverging 
low frequency response is killed as it comes entirely due to 
the gapless gauge fluctuations of the spin liquid. The res- 
onance due to the triplon ./V mode then becomes the most 
prominent feature in the (tt, tt) response. Further its fre- 
quency will soften as the doping is reduced. Second as the 
fermionic quasiparticles no longer have nodes at (tt, tt) 
they only weakly damp out this reonance. Finally the 
fermionic quasiparticles will still contribute some back- 
ground magnetic response which can now be usefully ad- 
dressed in a standard RPA calculation. Such calculations 
have been reported before in the literature, and give rise 
to incommensurate continuum scattering at frequencies 
below the resonance that appear to be consistent with 
experiment. 

Since the original discovery of the neutron resonance 
peak, there have been two more or less independent in- 
terpretations. One view is to describe it as a soft mode 
associated with the magnetism of the undoped Mott in- 
sulator. This view has the advantage that it provides 



a natural explanation of the softening of the resonance 
frequency with underdoping. The other view has been 
to simply regard it as a S — 1 collective mode of weakly 
correlated fermionic quasiparticles in the superconduct- 
ing state. The description given above unifies these two 
different interpretations. Indeed in the parent spin liq- 
uid the triplon N mode may be viewed as a particle-hole 
triplet exciton made out of spinons - rather than elec- 
trons. This mode appears naturally as a recognizable 
peak in the magnetic response upon approaching the AF 
state. Doping this spin liquid then leads to a supercon- 
ductor with gapless fermionic quasiparticles and a sharp 
gapped S = 1 triplon. 



VII. CONCLUSION 

In this paper we have revisited the issue of possible 
second order phase transitions out of the collinear Neel 
state into paramagnetic spin liquid states in two dimen- 
sional quantum antiferromagnets. The particular spin 
liquid we considered is a dRVB state which has gapless 
spin excitations. Correspondingly there are non-trivial 
power law corelations in the spin and other quantities. A 
useful description is provided in terms of gapless Dirac- 
like spinons that are coupled to a fluctuating U(l) gauge 
field. However there is possibly no true quasiparticle de- 
scription of the spectrum. Indeed this state is in a critical 
phase that is the two dimensional analog of the one di- 
mensional spin- 1/2 chain. In contrast to other simpler 
spin liquids which have a spin gap, a direct second order 
transition to the collinear Neel state appears to be possi- 
ble for such a two dimensional algebraic spin liquid. We 
developed in some detail a theory for such a transition. 
Magnetic long range order was obtained as a spin den- 
sity wave transition of spinons. We argued that gauge 
fluctuations convert the resulting magnetic state into a 
conventional one that is smoothly connected to the usual 
Neel state. Thus the spinons disappear from the spec- 
trum in the magnetic state. The theory for the transi- 
tion shares a number of similarities with the deconfined 
critical points studied recently. Most importantly, there 
are two diverging length/time scales as the transition is 
approached from the magnetic side. The shorter of the 
two scales is associated with the onset of magnetic order 
from a critical soup of spinons. The second longer scale 
is associated with confinement of the spinons. It is in 
the intermediate length/time scale regime (i.e. between 
the two diverging lengths) that the magnetic ordering is 
correctly described as a spin density wave formed out 
of spinons. This intermediate scale regime may also be 
characterized as a fractionalized antiferromagnet. 

We noted several implications of our results for theo- 
ries of the cuprates that regard them as doped spin liq- 
uids. First it allows us to develop a qualitative picture 
of the resonance mode seen in neutron experiments. Our 
picture unifies the existing descripions as a soft mode as- 
sociated with the magnetic ordering in the insulator and 
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as a triplet excition formed from a particle-hole pair of 
fcrmionic BCS quasiparticles. Indeed in our description 
the resonance is a soft mode of the magnetic ordering 
that is formed as a particle-hole triplet exciton of spinous. 
This picture is closest to that in RefS. 

We have shown how magnetism may be incorporated 
into the spin liquid based approach to the cuprates. Cen- 
tral to this is the description of magnetism as a spin den- 
sity wave ordering of spinons. Such a description has 
been explored before in a number of publications. As 
summarized in the first paragraph of this section, our 
work clarifies the range of validity of such a description. 
Indeed should experiments reveal clear signatures for a 
fcrmionic spinon description of the intermediate scale 
spin physics of the undoped cuprates then we could take 
that to be a signature of proximity to the quantum tran- 
sition to the dRVB algebraic spin liquid. 
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APPENDIX A: DIRAC ACTION 

Here we show that low energy effective action is de- 
scribed by a continuum theory of fermionc spinons with 
Dirac dispersion. As mentioned before, we need to 
expand the Hamiltonian close to nodes at (5,5) and 
(—§,§)• From now on (k x ,k y ) refer to deviation from 
(■|, J) °r (— § > f ) points. Here we explicitly derive the 
Hamiltonian near (^, ~). The other node is similar: 



H = 



2t k 4 



a 2 

+ 2iA fc_ 



2t fc. 



2iA fc_ 



9 2 



(Al) 



where k+ = k x + k y and fc_ = k x — k y . This Hamiltonian 
could be written in terms of two by two Pauli matrices 
which are defined as: 






1 


T y _ 


-i 


t z = 


1 


I— 1 







i 




-1 



Using this notation and also adding the contribution from 
(—5, § ) node we get the following form for low energy 
Hamiltonian: 

H = V c\ (2t k+r x + 2A k_T y - - t z N. — ) ci 
— g 2 



fe_L ,k- 



4 (2t k_T x + 2A k+r y + - t z N.^) c 2 

9 2 

(A2) 



Here c% and c 2 refer to (■§,§) and (—§,§) nodes, re- 
spectively. They are two component fermionic operators, 
each component representing one of the sites in the unit 
cell. The t 1 matrices operate in the space of these two 
components. Each component has a SU(2) spin index, 
where a matrices operates. Now assume t = A (i.e. ig- 
nore the velocity un-isotropy which is irrelevant in renor- 
maliztion group language) and we rename k + as k y and 
fc_ as k x . Then introduce the new fermionic operators: 



and subsequently: 



ipi = —ir x a 
tp2 = e 4 T r * c 2 



^1,2 = Wl,S rZ ) 



(A3) 
(A4) 



(A5) 



With these new variables, the Hamiltonian (|A2|) takes 
the following simple form: 



H = i>i{k x r x + k yT y + JN.^) fa 



(A6) 



+ Mk x r x + k yT y - jn.-) v> 2 



Then using fi matrices that operates in the space, made 
by presence of two different nodes, we can put but nodes 
contribution (ipi and ^2) as a single vector, ip. Then the 
Hamiltonian takes the form: 

H= $( k * rX + k v rV + U V Z #-|) ^ ( A7 ) 

This form in continuum limit and in real space leads to 
action given in 1121 



APPENDIX B: RANDOM PHASE 
CALCULATION 

We start with the partition function with an external 
magnetic field h that couples to the Neel vector: 



Z = J [Di/j][Di/)][Da][DN] e~ s (Bl) 

mixinq + S h (B2) 

S N = J d 2 xdT{^N) 2 + ^N 2 } (B3) 

S^,a = I (fxdTipi-i^idv+iea^ (B4) 



Smixing = J d 2 xdr iipa^i z ip.(XN) (B5) 
S h = [ d 2 xdTh.(ai{pafi z i} + bN) (B6) 



Here a and b are non-universal constants that depend on 
details of the microscopic physics. 
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Now by expanding the the terms which contain iip/j, z aip 
(one coming from S m ixing an d one from coupling to h 
field.) up to quadratic order, and performing the integral 
over ip and a fields we get the following form for the 
partition function: 




FIG. 8: Spin-Spin correlation in mean field 



Z = j[ DN ] e ~S^bf d ^drh.N J d 2 qduJ X^l\f l+X ^ 

(B7) 

Where x((j, cj) is the susceptibility of iipa^ip operator 
calculated using the action S^ ta . We can rewrite the 
x ( q ' u > \h + \N\ 2 up to quadratic order in N as an expo- 
nential. Putting this form in IB 71 we get the following for 
the partition function: 



lowing integral: 

(SiSj) = —5ij 



qeff 
N 



= Sn + / d 2 qduj{b h q .N- q — 



= I d 2 qdu>{ 



\Na 



2 

2 + h„.Nl 



\h + XN\ 2 } 



(B9) 



Now to get the quadratic action for h field we should 
integrate out N field, which is a simple Gaussian inte- 
gral now. This leads to the following value for effective 
susceptibility: 



Xeff 



a 2 X + b 2 Xn - 2Aa b \ Xn 
1 - A 2 x xn 



(BIO) 



You can see that as one expects, for A = 0, this reduces 
to some of susceptibilities for N and iip^, z ^ip fields. Also 
note that, this is valid only for the frequencies where 
\ 2 X Xn is not of order one {i.e. it is not valid for small 
frequencies). Now by analytic continuation to real fre- 
quencies and taking the imaginary part, we get the spec- 
tral functions plotted in section ITlI Al 



APPENDIX C: FEYNMAN DIAGRAMS 
1. Spin-Spin correlation in mean-field theory 

For spin-spin correlation we need to calculate the fol- 
lowing diagram which consists of two fermionic propaga- 
tors. The fermionic propagator in mean field is derived 
from the action [F2l 



i> - iXN.a 



(CI) 



This diagram, in the real space, corresponds to the fol- 



d 3 p d 3 q 
(2tt) 3 (2tt) 3 

-m. M(j+ 4-i\N.a){4-i\Nff)] 
((p + q) 2 + \ 2 N 2 )(q 2 + \ 2 N 2 ) 

(C2) 



(B8) After integrating over q this gives: 



i g + \N q \ 2 } 



(SiSj 



oc 



r 



p dp sin(pr) [4AJV+ 



- arctan(— — - 
p 2XN 



){4X 2 N 2 +p 2 )} 



(C3) 



The first term is proportional to 5(r). Since we are inter- 
ested in the case where r ^ 0, we can ignore that term. 
The second term gives: 



-2\Nr 



-(1 + 2AAV) 



(C4) 



We already saw that N goes to zero as g approaches g c 
like g c — g. So if we define correlation length as £ = ^Xn 
we get the equations IT31 and ITU 



2. Beyond mean field: e-expansion 

Now to do the e-expansion we need to study the ac- 
tion given in^| We can see that there are three different 
types of vertices present in the theory and eleven differ- 
ent one loop diagrams which cause field renormalizations, 
vertex corrections and mass renormalization. we use the 
following diagrammatic representation for the propaga- 
tors of the fields, in the theory: 



§ ab 



p 



P 
l 



(C5) 



Here we have used Feynman gauge for the gauge field 
propagator, and Euclidian metric. Now we start to cal- 
culat the one loop diagrams. 



14 



p - k 



J- 



p- k 

k 

(b) 



FIG. 9: fermion self energy 



c. Gauge field self energy 

Gauge field is only coupled to fermion field as well, 
so there is only one diagram generates gauge field self 
energy: 




a. Fermion self energy 

There are two, one loop dia grams which will generated 
fermions self energy. Diagram |9(a)| represent N field con- 
tribution to fermion self energy: 



£& i = (iA) 



d d k 



5 ij fi 



(2ir) d k 2 (k+p) 2 
3 A 2 , 



(C6) 



(4tt) 2 e 



Diagram 9(b) represents gauge field contribution to 



fermion self energ 



FIG. 11: Gauge Field Self Energy 
This diagram has the following contribution: 

d d k tr{Y WW+ f)) 



S5[p^-p^H = (-l)e 



(27r) d k 2 {k+pf 



16 e' 



p 2 [5^ v - 



(C9) 



Note that in this relation the term proportional to p^p 11 
dose not contribute to physical observables (S-matrix el- 
ements). This is guarantied by Ward identit}^. 



d d k u ft „ 

(27r) d7 k 2 (k+p) 21 
1 e 2 , 



(47T) 



(C7) 



ci. u((N) 2 ) 2 vertex correction 

There are two diagrams contributing to this vertex 
renormalization at one loop level. Note that up to one 



b. N Field Self Energy 

Since N field is coupled to fermion field only, there is 
just one diagram contributing to its self energy. 

P+k 

-;~0—- 



FIG. 10: N field self energy 



This has the following contribution: 

E^«F=(-l)(zA) J 



8 A : 
"(4tt) 2 e 



p 2 S ij 



(C8) 



Here the minus sign appears because of presence of a 
fermionic loop in the diagram^. 



k 



X 
1 



k 4 



(a) (b) 
FIG. 12: u renormalization 

loop order, it is enough to calculate ve rtex co rrections 
with zero external momentum. Diagram 12(a) gives: 

A%f(6 ij 5 kl +S ik S jl + S il S jk ) 

11? . . r d d k S oq S pt 

= (s l 3 6 op + S^d™ + 5 ip 5 JO ) / : — 

6 { 'J (2tt) 2 (fc 2 +r) 2 



11 u z 



(S qt 5 kl + S qk S tl + 8 ql 8 tk ) 
{5 t3 8 kl + 5 lk 5 jl + 6 il S jk ) 



(CIO) 
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The contribution of diagram 12(b) is similarly calculated: 



A%(S ij 8 kl + 5 ik 6 jl +5 il 6 jk ) 

= -^-^V ?S^'* vW) (on) 



96 X 4 
(4^7 



(5 ij S kl + S ik S jl + S il 6 jk ) 



Again one minus sign appears because of the fermionic 
loopi^. 



e. N-ip vertex correction 

Again, there are two different diagrams associated with 
this correction at one-loop level: 





(a) 



(b) 



FIG. 13: A renormalization 



As before we set the external momentums to zero. Then 

is: 



the contribution of diagram 13 (a 
iA x N a k = (iX) 3 



d d k ft i S ij a k ^ ft 
(2^F 7 k 2 +r ° fc 2 
2 A 3 , 



(C12) 



and diagram 13(b) 



(4tt) 2 e 
gives: 



iA* =iX e z 



d d k ft „ §" v u ft 



(2ir) d k 2 7 k 2 +r' k 2 



I 7,2 



6 A 



(C13) 



(4tt) 2 e 



/. a-ip vertex correction 

Similar to lC 2 el there are two diagrams: 
The first diagram 1 14(a) I gives: 



d d k u ft 7°^ ft „ 



(27r) d7 k 2 k 2 k 2 

e 3 „ 



(C14) 



(An) 2 e 



■1 





(a) 



(b) 



FIG. 14: e renormalization 



and 14(b) gives: 

r =e(zA) 2 

3 e 3 
~(4tt) 2 e 



d d k t ft j a S^ ft 



(2n) d k 2 k 2 + r k 2 



(C15) 



'7 



Note that EH is minus the contribution of IC7I and lC15l 
has the minus contribution of IC6I This, in fact, should 
be the case to keep gauge invariance of the theory. 



g. Mass renormalization 

The last one loop diagram, we consider, generates mass 
renormalization: 

k x / 1 



A r S' j = -(S' j 5 kt + 5 ik S jl + S il 6 jk ) 
6 



d d k d 
(27tFF" 



Ik 



5 u 
(4tt) 2 37 



6 lJ 



(C16) 



3. Renormalization conditions 



Now with these in hand, we can proceed using minimal 
subtraction scheme. Introducing a mass scale m oc ^/r we 
can write the following set of renormalization conditions: 



Z N (p 2 +r] 



^Z~ a + (A e N + Al) 
+ (rA r ££ p 2 ) 
m-iXoZtpy^ + (A X N + A„) 
P 2 (Z a ~ 
m- e u Z 2 N + (A« + AJ) 



finite 0(2 loops) 
finite 0(2 loops) 
finite 0(2 loops) 
finite 0(2 loops) 
finite 0(2 loops) 
finite 0(2 loops) 



16 



Here eo, Ao and uq are bare coupling constants and so 
does not flow with mass scale. These relations give field 
renormalization coefficients directly, since the divergence 
part of self energy diagrams should cancel out with these 
field renormalization coefficients. 



Zn 

Z a 



N 



Ea 



(C17) 
(C18) 
(C19) 



Putting them back in renocmalization condition equa- 
tions, and letting the mass scale fiowiS, we get the equa- 
tions given in section Ivl 



4. Renormalization of bilinear operatores 

Here we introduce a new term in the Lagrangian with 
the general form: 



(C20) 



where O is the combination of /i and a matrices which 
are given in table [i] There are two new diagrams corre- 
sponding to this new vertex ( 15(a)|15(b) I Diagram 15(a) 





This gives the scaling dimension of v (1 + S v ). Now to 
get the scaling dimension of O (Ao), note: 



A 



o 



D-(l + 6 v ) 



(C24) 



Using the fixed point values for A and e, we get the scaling 
dimensions mentioned in table [i] 



5. Velocity anisotropy 

In this section we assume a small velocity anisotropy 
and treat it as a perturbation to our QED3 theory. Fol- 
lowing the notation used in Ref. l30l this anisotropy is 
presented by: 



K„ 



(C25) 




(a) 



(b) 



FIG. 16: Velocity anisotropy vertexes 



(a) 



(b) 



gives: 



FIG. 15: bilinear operators scaling dimension 



„i,,z ft rn ft ^ ,,z j 



k 2 k 2 



A o 0, A = o l v z Op J z a l 



(4tt) 2 e 



and diagram 15(b) 



A? C =e 2 



gives: 
d d k 

(2n) d 1 k 2 ~ k 2 k 2 



' 1.9 ^ 1.?. I."?. I 



G 



(C21) 



(C22) 



(4tt) 2 e 



O 



Whit these in hand we can get the scaling dimension of 



A„ = 1 + S v = [1 + (Ag + )A 2 + (A a a + 

A 2 e 2 
- [1 + (-2A - 3)-— j + 5— -nrl» 



(4tt) 2 (4tt) s 



(C23) 



Here, S is the small perturbation parameter (we are es- 
sential interested in its behavior under renormalization) 
and 7^ = 7 K £a; )(L1 — j v S y ^. There are two vertexes as- 
sociated with this perturbation. One is correction in 
fermionic kinetic energy presented in figure 16(a)| and 
correction to the a — ip vertex presented in figure |16(b) " 



-5n z }t 



(C26) 
(C27) 



So there are there three different type of one-loop 
diagrams, in addition to field renormalization factors 
calculated before contributing renormalization of the 
anisotropy term in fermionic kinetic energy {i.e. S). Di- 
agram 17(a) gives: 



p- k 



p - k p- k 



k 

(a) 



k 
00 



k 

(c) 



FIG. 17: Velocity anisotropy one-loop renormalization 
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-6 E$tf=-5 {i\y 



(2n) 



ji{p-k) 2 + r 



S A 2 
(4tt) 2 e 



(C28) 



Now with these in and also the field renormalziation co- 
efficients calculated perviously we can get the RG flow 
for 5: 



Diagram 17(b) similarly gives: 



-5W 



-5e< 



3(4tt) 2 3e 



d d k „ 1 - 1 6^ _ 



(C29) 



ft = (5 (2 + 



3(4^) 2 



(4tt) 5 



(C31) 



Diagram 17(c) and the similar one with the correction on 
the right vertex give: 



-8 A a s ^ = S e 2 / — — r 



(27r) d ' ft (p-k) 2 



S 2e 2 
(4tt) 2 e 



(C30) 



This proves that the small velocity anisotropy is irrele- 
vant. 
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